Metastable Vortex Lattice Phases in Superconducting MgB2 
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The vortex lattice (VL) symmetry and orientation in clean type-II superconductors depends 
sensitively on the host material anisotropy, vortex density and temperature, frequently leading to 
rich phase diagrams. Typically, a well-ordered VL is taken to imply a ground state configuration for 
the vortex-vortex interaction. Using neutron scattering we studied the VL in MgB2 for a number 
of field-temperature histories, discovering an unprecedented degree of metastability in connection 
with a known, second-order rotation transition. This allows, for the first time, structural studies 
of a well-ordered, non-equilibrium VL. While the mechanism responsible for the longevity of the 
metastable states is not resolved, we speculate it is due to a jamming of VL domains, preventing a 
rotation to the ground state orientation. 

PACS numbers: 74.25. Ha, 74.70.Ad, 74.25.Wx, 61.05.fg 



Metastable phases of matter are well-known, with fa- 
mous examples including supercooling and superheating 
of liquids and diamond which is one of the many al- 
lotropes of carbon. Metastability is almost exclusively 
observed in connection with first-order transitions, and 
is often found in frustrated systems where the energy dif- 
ference between the states is small. The structure of the 
vortex lattice (VL) in typc-II superconductors is known 
to be highly sensitive to changes in external parameters 
such as temperature and magnetic field and can therefore 
naively be expected to display metastability, for example, 
in connection with first-order transitions such as the VL 
melting [l[ or the reorientation transition of the rhom- 
bic VL found in most superconductors with a four-fold 
in-plane anisotropy However, except for phase coexis- 
tence in the vicinity of reorientation transitions , no 
well-ordered metastable VL phases have been observed 
to date. 

We present extensive small-angle neutron scattering 
(SANS) studies of the VL in MgBz with H || c which 
is found to exhibit an unprecedented degree of metasta- 
bility. Previous studies of this compound revealed a con- 
tinuous, field-driven 30° VL rotation bounded by two 
second order transitions. The rotation is understood to 
arise from a competition between six-fold Fermi surface 
anisotropics with opposite signs on the tt- and cr-bands, 
coupled with the suppression of the 7r-band supercon- 
ductivity by a modest magnetic field 0, [1]. The most 
striking new results are the existence of highly ordered 
non-equilibrium VL configurations and the observation of 
hysteresis in connection with a second order phase tran- 
sition which we speculate is the result of a jamming-like 
scenario. As a result, the ground state VL phase diagram 
is found to differ significantly from our earlier report as 
well as theoretical predictions 0, S] ■ 



Small-angle neutron scattering (SANS) experiments 
were performed at the D22 instrument at Institut Laue- 
Langevin, the HFIR GP-SANS beam line at Oak Ridge 
National Laboratory and NG3 at the NIST Center for 
Neutron Research, using a standard configuration with 
the applied magnetic field parallel to the incoming neu- 
tron beam [91] . To achieve diffraction the sample and the 
magnet were rotated and/or tilted together in order to 
satisfy the Bragg condition for the VL planes. Complete 
VL diffraction patterns were obtained by adding the scat- 
tering, recorded by a 2D position sensitive detector, from 
several angles. To resolve closely located VL Bragg re- 
flections, a very tight coUimation of the neutron beam 
was used (D22: 0.06°, GP-SANS/NG3: 0.08° FWHM). 

Several different single crystal samples were studied 
providing consistent results, with only slight variations 
in the VL transition fields/temperatures. The crystals, 
with masses in the range 100 - 200 fig, were grown using 
isotopically enriched ^^B to reduce neutron absorption 
and were comparable to those used in previous SANS 
experiments by our group 0, EB, 11 1 ■ 

Measurements were performed at temperatures be- 
tween T = 2 and 28 K, and in magnetic fields up to 
iJ = 1.9 T applied parallel to the crystalline c axis. Sev- 
eral different field-temperature histories were employed: 
Constant field cooled (FC) from a temperature above Tc 
and field ramped (FR) where H was changed at a con- 
stant temperature. Other measurements were performed 
following a damped, small-amplitude field oscillation ap- 
plied after either a FC or FR procedure. 

Fig.[l]shows VL diffraction patterns, obtained at three 
different places in the mixed state phase diagram of 
MgB2, which summarize the main result of this report. 
Panels (a) and (b) were obtained following a field cool- 
ing (FC) procedure from a temperature T > Tc to 18 K 
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FIG. 1. SANS diffraction patterns showing the different triangular VL phases in MgB2 with H || c. The orientation of the 
crystalline a axis is shown in panel (a) and the a* axis in panel (e). The three columns show the VL configurations at 18 K 
and 0.5 T (a and d), 2 K and 0.5 T (b and c), and 2 K and 1.1 T (e and f). In cases where not all VL reflections were rocked 
through the Bragg condition, symmetry equivalent positions, obtained by mirroring through and/or rotating around q ~ 0, are 
indicated by open circles. In the high symmetry F- and I-phases a single-domain VL is observed aligned with Bragg reflections 
along respectively a and a*. In the intermediate L-phase two VL domains orientations, indicated by 2 x 6 = 12 peaks, are 
observed rotated away from a by an angle Lp in opposite directions as shown in panel (f). For each field/temperature the ground 
state VL is obtained following a damped field oscillation with an initial amplitude of 25 mT. 



and then 2 K in a constant field of 0.5 T, and both show 
VLs oriented with Bragg reflections along the crystalline 
a axis. In both cases the VL is well-ordered, indicated 
by very sharp and well-defined Bragg peaks and narrow 
rocking curves shown in the Supplemental Material. Tra- 
ditionally such high quality diffraction patterns are taken 
as an indication that the VL is an equilibrium configura- 
tion for the vortex- vortex interaction However, ap- 
plying a small-amplitude field oscillation at 2 K causes a 
transition to a VL phase which is rotated away from the 
a axis by an angle ip, as shown in panel (c). Here a super- 
position of scattering from two degenerate domain orien- 
tations is observed, corresponding to clockwise and anti- 
clockwise rotations of the VL. The widths of the peaks 
(radial, azimuthal, longitudinal/rocking curve) stay effec- 
tively unchanged. The rotated VL remains robust upon 
heating to 18 K as shown in (d), but a subsequent field 
oscillation causes a transition to the a axis aligned phase 
as in panel (a). Additional measurements showed that in 
the absence of a field oscillation the rotated VL phase at 
0.5 T (d) persists up to the highest measurable temper- 
ature of 28 K. 

These results show that the ordered VL in MgB2 ex- 
hibits mctastability as it is heated or cooled across the 



equilibrium phase transition which must exist somewhere 
between 2 and 18 K. Furthermore, it is necessary to in- 
duce vortex motion by a small- amplitude field oscillation 
in order to drive the VL to the ground state (equilib- 
rium) configuration. A similar effect is observed at 1.1 T 
where the field cooled VL is aligned along the a* axis as 
shown in panel (e). After a field oscillation, the VL enters 
the rotated phase (f). In summary, this system presents, 
for the first time, the opportunity to perform structural 
studies of a well-ordered, non-ground state VL. Due to 
the similarity with the tilted hexatic phases in liquid crys- 
tals, we adopt the same naming scheme, denoting the VL 
phases as F, I and L as indicated in Fig. [T] [l2- 14|. 

Fig. [2][a) shows a comparison between the VL orien- 
tation angle ip, obtained at 2 K following respectively a 
FC and a field ramping (FR) procedure where the mag- 
netic field is changed at the measurement temperature. 
For fields > 0.2 T the two procedures result in differ- 
ent VL phases. Furthermore, the ground state VL at 
this temperature never completes the rotation transition 
but remains in the L-phase, stabilizing at « 24° above 
1.2 T. We note that a FR or a FC followed by a field 
oscillation both achieve the same result. Fig. [2Kb) shows 
the field dependence of the FR VL orientation obtained 
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FIG. 2. (Color online) VL orientation angle as a function of 
applied field. Panel (a) shows the the results of measurements 
at 2 K following a FC procedure and a zero field cooling fol- 
lowed by a FR. The arrow indicates the FL-transition for the 
ground state VL. Panel (b) shows the FR VL orientation at 
temperatures ranging from 2 to 20 K. All lines are guides to 
the eye. 

at temperatures in the range 2 - 20 K. Only for temper- 
atures of 10 K and above does the ground state VL reach 
the I-phase {^p = 30°). The ground state VL configura- 
tion is summarized in the phase diagram in Fig. [31 Due 
to decreasing intensity, it was not possible to determine 
the phase boundaries above 1.8 T and 20 K. Nonetheless, 
the constant for iif > 1.3 T at both 2 and 5 K seen in 
Fig. [2] strongly suggests that Ll-transition is essentially 
vertical at low temperatures as indicated by the extrap- 
olation in the phase diagram. Likewise, the fact that the 
L-phasc can be observed following a FC implies that it 
extends up to Hc2 within a narrow field range. 

The ground state VL phase diagram differs markedly 
from a naive expectation that the effect of the 7r-band 
anisotropy, believed to be responsible for the F- and L- 
phases, should be largest at low fields (little suppres- 
sion of TT-band superconductivity) and high tempera- 
tures (increased thermal mixing of bands) 0, Hi- Con- 
sequently, one would expect the I-phase to occupy the 
high- field/low-temperature part of the phase diagram in 
striking contrast to Fig. [3l In particular, the reappear- 
ance of the L-phase as one follows Hc2 to lower temper- 
atures is unexpected and does not agree with theoretical 
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FIG. 3. Ground state VL phase diagram for MgB2 with 
H II c. The FL-transition is indicated by squares and the 
Ll-transition by circles. The lines are guides to the eye with 
dashes indicating extrapolations to higher temperatures and 
fields. Diamonds indicate measured values of Hc2- 

predictions [Ij. 

We now return to the metastable VL phases shown in 
Fig. [IJb), (d) and (e), noting that the existence of such 
states in connection with a second order phase transi- 
tion is unexpected as shall be discussed in the follow- 
ing. It should be stressed that the present case is not 
one of a first order transition broadened by sample inho- 
mogeneities as for example observed in connection with 
the order-disorder transition associated with the peak ef- 
fect [l5j. As pointed out by Zhitomirsky and Dao, the 
continuous rotation of the VL in MgB2 with second order 
FL- and Ll-transitions shows that the VL free energy as 
a function of the orientation angle must be given by 

5F = Kq cos 6ip + Ki2 cos 12(p, (1) 

with Ki2 > 0. The same free energy expression 
has been used to model the reorientation transitions 
observed between tilted hexatic phases in liquid crys- 
tals [iHl, where we note no metastable phases are ob- 
served. Fig. lU^a) shows the evolution of 5F as Kq is 
increased from < — 4ifi2 to > 4i^i2. In Fig. |4ljb) we 
consider the free energy curves corresponding to the cy- 
cle of VL phases in Fig. [IJa - d), which shows that 
the metastable F-phase (b) corresponds to an unstable 
equilibrium [d{5F)/dip = 0; (f{SF)/d(p^ < 0] while the 
metastable L-phase (d) represents a true non-equilibrium 
situation [d{SF)/d(p ^ 0]. This implies that the VL 
metastability can not be understood from the single do- 
main free energy thus raising the important question: 
Which mechanism is responsible for the longevity of the 
metastable states, preventing them from immediately ro- 
tating to the ground state? 

A straightforward answer to the above question is that 
the metastable VLs are held in place by vortex pin- 
ning. However, in general metastability is a hallmark 
of very clean systems, in obvious contrast to such a see- 
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FIG. 4. (Color online) Free energy curves for the VL phases in 
MgB2. Panel (a) shows the progression of free energy curves 
calculated using Eqn. ([1} with K12 = 1 and increasing Kq. 
The corresponding ground state VL phase is indicated by each 
curve. The FL- and Ll-transitions occur at Kq = ±4. Panel 
(b) shows free energy curves corresponding to the sequence 
of diffractions patterns at 0.5 T shown in Fig.[T] Solid circles 
represent ground state VL configurations while open circles 
represent metastable phases. 



nario. More specifically this also seems unlikely given 
the very low depinning critical current in clean MgB2 
single crystals similar to the ones used in the present 
work (< 2 X 10^ A/cm^ at 5 K and 0.5 T and decreas- 
ing with higher T and H) [13]. Moreover, at 0.5 T the 
irreversibility line for H || c is in the range 8 - 14 K [Tt} . 
much below the highest measurable temperature of 28 K, 
where we still observe a metastable L-phasc. Finally, the 
sharp VL Bragg peaks observed in the FC case does not 
indicate any significant disorder due to vortex pinning in 
contrast to reports for e.g. YNi2B2C 0-0] ■ 

We speculate that the mechanism which sustains 
the metastable configurations is a jammin g o f counter- 
rotated VL domains against one another [18|. In this 
scenario the VL domains acts as granular entities, pre- 
venting each other from rotating to the ground state ori- 
entation. If this is the case a close examination of the 
metastable F- and I-phases should find an onset of do- 
main formation, leading to an azimuthal broadening of 
the VL peaks which would be most visible at higher fields 
as the resolution improves with H . Such an azimuthal 



broadening is indeed observed at 2 K and 1.1 T, increas- 
ing from a resolution limited ~ 4° in the ground state 
L-phase to ^ 7.5° in the metastable I-phase, as shown 
in Fig. S3 and discussed in detail in the Supplemental 
Material. No broadening is expected or observed in the 
metastable L-phase. Further studies are needed to fully 
explore this jamming hypothesis. A key question in this 
regard will be to investigate the robustness of the VL 
domain boundaries necessary to prevent rotation. Theo- 
retically, recent work by Deutsch and Shapiro considered 
domain boundaries in a square VL [l^. Such studies 
should be extended to investigate how this affects the 
VL configuration. 

Before concluding, we note that SANS studies of the 
VL in UPta used a temperature quenching technique 
to infer information about the superconducting A-phase 
which exist just below Hc2 in this material |20|. In this 
case a rotated VL identical to the L-phase was inter- 
preted as originating from a particular order parameter 
symmetry. Our results on MgB2, which does not have 
the exotic pairing symmetry of UPta, are qualitatively 
similar, suggesting that this interpretation may have to 
be reconsidered. 

In summary, we discovered a high degree of metastabil- 
ity in the MgB2 VL, which allows measurements in non- 
equilibrium configurations. The ground-state VL phase 
diagram was determined and found to differ substantially 
from theoretical predictions. We speculate that a jam- 
ming scenario is responsible for the VL metastability. 
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Here we present an analysis of the MgB2 vortex lat- 
tice (VL) Bragg peak widths measured by small-angle 
neutron scattering (SANS) with H || c. The measure- 
ments were performed following a variety of magnetic 
ficld/tcmpcraturc histories to explore both ground state 
and metastable VL configurations. Measurements of the 
ground state phases were performed following a magnetic 
field ramp at constant temperature with both increasing 
(FR up) and decreasing (FR down) fields, or following a 
small-amplitude (25 mT) damped field oscillation (FO) 
after a temperature change. Metastable VL phases were 
measured when a temperature change across one of the 
equilibrium phase transition was not followed by a FR or 
FO. 

Three widths can be measured, as indicated in Fig. [1] 
The radial {Sq-s) and azimuthal ((^v^a) widths are related 
to the positional and orientational order of the VL, while 
the longitudinal curve width (iJwl) is a measure of the 
straightness of the vortices. The first two widths were ob- 
tained by a two-dimensional Gaussian fit to the VL peaks 
recorded on the detector (Fig. 1 in main paper). The 
longitudinal width was obtained from so-called rocking 
curves, which is the intensity variation of the VL Bragg 
peaks as they are rotated through the Ewald sphere. The 
rocking curves were corrected for the Lorentz factor (an- 
gle at which peaks move through the Ewald sphere) and 
fitted by a Voigt function. Further details concerning 
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SANS and how the rocking curves are measured can be 
found in Ref. ? . 

The measured widths will be compared to the ex- 
perimental resolution obtained by the following expres- 
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Here k = 27r/A„ is the neutron wavevector magnitude 
and A„ is the wavelength. The VL scattering vector 
(for a hexagonal symmetry) is given hy q = 1.075 x 
27r(_B/0o)^/^ where B is the magnetic field and 4>q — 
h/2e — 2067 Tnm^ is the flux quantum. There are two 
contributions to the resolution. The first comes from the 
wavelength spread AA„/A„ = 10% and is given by 
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FIG. 1. Schematic showing the three different VL widths 
measurable by SANS. 



where the numerical factor is due to the conversion from 
variance to FWHM. The second contribution (A(20)bd) 
is from the beam divergence, which was determined from 
the size of the undifFractcd beam spot on the detector. 
Since k ^ q the longitudinal resolution is substantially 
higher than the in-plane resolution. 

Fig. [2] shows the field dependence of the radial width 
at 2 K. This is found to increase with field, but re- 
mains resolution limited throughout the range of mea- 
surements. In particular there is no measurable differ- 
ence between ground state and metastable VLs, showing 
that within the measurement precision the variation of 
the lattice plane spacing at all fields is the same for both 
phases. The calculated experimental resolution is possi- 
bly slightly overestimated. The measured widths at high 
field, which fall substantially below the resolution, can 
be attributed to the rapidly decreasing scattering inten- 
sity and resulting poor fits which yield unphysical results. 
The radial width can, in principle, be used to estimate 
the size of the VL domains in the direction perpendic- 
ular to the field. However, since this width is resolu- 
tion limited {SqB. ~ Agji) it only provides a lower bound 
Cr > 2/SqH « fOO - 200 nm. 

We now turn to the field dependence of the azimuthal 
width shown in Fig. [3l Here all ground state VLs 
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FIG. 2. Field dependence of the radial width ((5r) of the 
VL Bragg peaks at 2 K for different field/temperature histo- 
ries: Increasing field ramp (FR up), decreasing field ramp (FR 
down), field oscillation (FO) and cooling/heating in constant 
field (Metastable). The solid line is the experimental resolu- 
tion with the discontinuities corresponding to changes of the 
neutron wavelength (A„ = 10, 7 and 5 nm with increasing 
field). 
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FIG. 3. Field dependence of the azimuthal width ((5<^a) of 
the VL Bragg peaks at 2 K. The solid line is the experimental 
resolution with the discontinuities corresponding to changes 
of the neutron wavelength. 



are again found to be resolution limited, whereas the 
mestable VL at 1 T is found to be significantly broad- 
ened. Similarly the metastable VL at 0.5 T is broader 
than the ground state but in this case this is within 
the statistical fluctuations of the data. For the ground 
state we find a lower bound of the VL domain size of 
Ca > 2/((7(5(y5A) ~ 150 — 300 nm in good agreement with 
the estimate or above. We interpret the broadening of 
the metastable VL at 1.1 T as being due to a small rota- 
tion (onset of transition to L-phase) which is insufficient 
to clearly separate the peak into two. 

The field dependence of the VL rocking curve width 
is given in Fig. 21 This show a conventional monotonic 



decrease with increasing field, indicative of a straight- 
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FIG. 4. Field dependence of the rocking curve width (&i>Ji-C) of 
the VL Bragg peaks at 2 K. The solid line is the experimental 
resolution with the discontinuities corresponding to changes 
of the neutron wavelength. 
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FIG. 5. Temperature dependence of the radial width of the 
Bragg peaks at 0.5 T in the ground state (FO) and metastable 
VL phases. The solid line is the calculated experimental res- 
olution for A,i = 7 nm. 



ening of the vortices as they are pushed closer together 
and their mutual repulsion become stronger. At roughly 
1 T the rocking curve width becomes resolution limited. 
An estimate of the length of straight vortex segments is 
given by > 2/(5 Jwl) which varies from 4 to 6 /im. 
No difference between metastable and ground state VL 
configurations is observed. 

Figs. [5] to [7] show the temperature dependence at 0.5 T 
of the radial, azimuthal and rocking curve widths respec- 
tively. In all cases no temperature dependence is ob- 
served within the measurement precision. Similarly, no 
differences between metastable and ground state VLs are 
found. In the case of the radial width we again find this 
to be resolution limited, with fits at high temperature 
which fall below the experimental resolution. The low 
intensity at high temperature also leads to very large er- 
rors on the rocking curve widths as seen in Fig. [71 
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FIG. 6. Temperature dependence of the azimuthal width of 
the VL Bragg peaks at 0.5 T. The sohd hne is the experimen- 
tal resolution. 
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FIG. 7. Temperature dependence of the rocking curve width 
of the VL Bragg peaks at 0.5 T. The solid line is the experi- 
mental resolution. 



